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ABSTRACT 

Context. Solar spots appear to decay linearly proportional to their size. The decay rate of solar spots is directly related to magnetic 
diffusivity, which itself is a key quantity for the length of a magnetic-activity cycle. Is a linear spot decay also seen on other stars, and 
is this in agreement with the large range of solar and stellar activity cycle lengths? 

Aims. We investigate the evolution of starspots on the rapidly-rotating (Pj-ot ~ 24 d) KO giant XX Tri, using consecutive time-series 
Doppler images. Our aim is to obtain a well-sampled movie of the stellar surface over many years, and thereby detect and quantify a 
starspot decay law for further comparison with the Sun. 

Methods. We obtained continuous high-resolution and phase-resolved spectroscopy with the 1.2-m robotic STELLA telescope on 
Tenerife over six years, and these observations are ongoing. For each observing season, we obtained between 5 to 7 independent 
Doppler images, one per stellar rotation, making up a total of 36 maps. All images were reconstructed with our line-profile inversion 
code iMap. A wavelet analysis was implemented for denoising the line profiles. To quantify starspot area decay and growth, we match 
the observed images with simplified spot models based on a Monte Carlo approach. 

Results. It is shown that the surface of XX Tri is covered with large high-latitude and even polar spots and with occasional small 
equatorial spots. Just over the course of six years, we see a systematically changing spot distribution with various timescales 
and morphology, such as spot fragmentation and spot merging as well as spot decay and formation. An average linear decay of 
D - -0.022 + 0.002 SH/day is inferred. We found evidence of an active longitude in phase toward the (unseen) companion star. 
Furthermore, we detect a weak solar-like differential rotation with a surface shear of a = 0.016 + 0.003. From the decay rate, we 
determine a turbulent diffusivity of rjx - (6.3 + 0.5) x 10^"^ cm^/s and predict a magnetic activity cycle of ^ 26 + 6 years. Finally, we 
present a short movie of the spatially resolved surface of XX Tri. 
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1. Introduction 

Sunspots enable us to trace differential rotation and meridional 
circulation on the solar surface to extreme precision as well as 
to measure in detail a surface magnetic field in full three di¬ 
mensions (e.g.. lSolanyi2003l:lMoradi et al.ll201Ql) . Sunspots are 
understood as the emergence of magn etic flux tubes originating 
from a dynamo p rocess in the interior (lBabcocklll96ll: iLeightonI 
I^Ml989l) . in mean-field terms referred to as the so-called 
crQ-dynamo. As sunspots host strong magnetic fields, their study 
may provide indirect information about the internal dynamo ac- 
tivity, reminding us that the latter is still very controversial (e. 


^_ . erv c , 

iKanvla et al.ll2()rill2(nHTahhari et al.ll2fn4:lYadav et aTII^ 

Even the decay of sunspots, and thus the decay of 
the surface magnetic flux, is still not frUv understood 
(iRudiger & Kitchatinovll200Ql; iMartmez Pilletll20Q^ . Analyzing 
a subset of th e Greenwich Photoheliographic Results (GPR), 
iBumbal (Il963l) proposed a linear decay law of form dA/dt = D 
(where A is the area of the sunspot usually given in units of 
millionths of the solar hemisphere, 1 MSH = 3.05 Mm^) for re¬ 
current sunspots (those with two or more disk passages) with 
a mean value for D of -4.2 M SH/day and an exponential 
law for nonrecurrent spot groups. iMoreno-Insertis & VazquezI 


^ Based on data obtained with the STELLA robotic telescopes in 
Tenerife, an AIP facility jointly operated with I AC. 


(Il988l) analyzed the GPR data from 1874-1939 and obtained a 
parabolic decay law. Using the more accurate Debrecen data, 
iPetrovav & van Driel-Gesztelyil(ll997[) obtained re s ults su pport- 


mg a parabolic decay law. 
lyzed sunspot decay rates 


Martinez Pillet et all (1 19931) ana- 
using the GPR data from 1874- 


1976 and found a lognormal distribution as well as some evi¬ 
dence for weak nonlinearities in the decay process of isolated 
spots. A linear area de cay law has been studied from a the¬ 
oretical point of view ^okhale & ZwaanI [197^ Mever et al.l 
11974 iKrause & RudigeJ[l975l) . iMever et ^ (11974 ) obtained a 


constant area decrease rate through a process of diffusion of 
magnetic field over the entire area of the spot. This diffusion 
model predicts a linear area (and magnetic flux) decay, imply- 
ing that dAjdt is propo rtional to the turbulent diffusivity tjt- 
iKrause & Rudiged (119751) proposed a similar model based on tur¬ 
bulence, where the turbulent diffusion was related to the flux de¬ 
cay by d<f>ldt oc y/T. lRudiger & KitchatinovI (120001) used a mean- 
field formulation of diffusivity quenching and produced quasi- 
linear decays for both the spot area as well as the magnetic flux. 

Since the mo dern rediscovery of starspots (Eal 119721: 
IStrassmei^ l2009l and the many references therein), the study 
of spots on other stars introduced an extra line of research to 
better understand the Sun. Despite this, it is generally not possi¬ 
ble to reso lve stellar surfaces direct ly, with a few recent excep¬ 
tions (e.g. jKloppenborg et al.ll201Qh . clever mathematical meth- 
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Fig. 1: STELLA spectroscopy of XX Tri. a) A section of a representative single spectrum showing a subset of the spectral lines used 
for the Doppler imaging inversions, b) The full time-series radial velocities from a cross-correlation analysis. It shows the rotational 
phase sampling for the STELLA spectra. Each dot is from a full echelle spectrum, covering the wavelength range 390-880 nm. c) 
A section of b) referring to the observational season 2008/09. The dotted line represents a radial-velocity fit. 


ods and observing technique s were introduced to resolve stel¬ 
lar surfaces indirec tly (e.g., IVogt & PenrodI Il98"^ iRicd 1200^ 
ICarroll et"^ l2012h . These techniques, commonly referred to 
as Doppler imaging, became the most advanced tool for the 
study of starspots. Unfortunately, the needed time series of well- 
sampled, high-resolution spectra are difficult to obtain, even for 
a single Doppler image, given that intrinsic stellar surface vari¬ 
ations confine the data gathering process to a single or possibly 
two consecutive stellar rotations. For a star like XX Tri with a 
rotation period of 24 days, this is obviously not a simple task. 
Sampling the evolution of starspots requires repeated visits of 
a target, or preferably even continuous decade-long time series 
of spectra. In the present paper, we present for the first time, 
such long-term, highly-sampled, phase-resolved spectroscopic 
data for Doppler i maging, made possible by the use of STELLA 
robotic telescopes (IStrassmeier et al.|[2010bh . 


There have been several other attempts to monitor active 
stars by means of Doppler imaging. Eor example, for the rapidly 
rotating single late-type giant EK Comae, a series of snap¬ 
shots co nsisting of 25 Doppler maps between 1993-2003 ex¬ 
ists, see iKorhonen et al.l (l2007l). The RS CVn binary II Peg 
was monitored bv iBerdvugina et"^ (ll998L[T9^ between 1992- 


s. The se maps were par- 
20111) . including six new 


1999, resulting in 1 5 Doppler map s 

tially reanalyzed bv iLindborg et alT (f _ , _^_ 

maps. iHackman et al.l (l2012h added 12 more Doppler maps us¬ 
ing the same spectrograph (SOEIN at the Nordic Optical Tele¬ 
scope), resulting in a total of 28 Doppler maps between 1994- 
2010. Another RS CV n binary, IM Peg, was monitored by 
iBerdvugina et al.l ( 2000l) who obtained eight Doppler maps be¬ 
tween 1996-1999. iMarsden et al.l (l2007l) obtained a movie from 
31 Doppler maps for IM Peg during a monitoring campaign be¬ 
tween 2004-2007. Another starspot movie was published ear- 
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Her for HR 1099 (V711 Tau) from 37 Doppler maps taken 
in 199 6 (IStrassmeier & Bartua l200Qh . Table 2 in IStrassmei^ 
(l2009l) lists all Doppler-imaging attempts and their general re¬ 
sults until 2008. 

Our target is the spotted, red giant XX Tri (HD 12545), a 
member of the RS CVn class of magnetically active compo¬ 
nents of close binaries. This red giant is a synchronized SBl- 
type system with a period of Prot ~ ^orb ~ 24 days. We have ob¬ 
tained 667 usable spectra between July 2006 and April 2012, 
which cover 36 rotational periods. This star is famous for its de- 
tected superspot w ith a linear extension of 12x20 solar radii 
(IStrassmeiej 1 19991) . In the mid-1980’s XX Tri became one of 
the most attractive targets for phot ometrists because of its un¬ 
usually large lightcurve variations. iBidelmanI (Il985l) observed 
strong Ca ii H&K emission lines in the spectrum of XX Tri and 
suggested that it might also be a photometric variable. Photo¬ 
metric variability of RS CVn stars is generally attributed to the 
presence of large, cool starspots moving in and out of view as 
the star rotates. The presence of an extremely active chromo - 
sphere on XX Tri was confirmed by IStrassmeier et al.l (ll99Qh . 
showing that the Ca ii H&K emission intensity was two to three 
times that of the local continuum. The first published photomet¬ 
ric observations d ated back to 1986-87 and showed a fairly scat¬ 
tered light curve (iHooten & Halil [19^ . Three year s later, a re¬ 
marka bly high V amplitude of 0^6 was observed bv iNoltheniusI 
(Il99ll) . Multicolor photometry, obtained in early 1991, show¬ 
ing a large amplitude of 0?5 in V and 0?12 in V - and al¬ 
lowed IStrassmeier & Qlahl (Il992h to derive a precise tempera¬ 
ture difference between spot and photosphere of 1100 + 35 K, 
which suggests a spot coverage of approximately 20 % of the 
entire stellar surface. Furthermore, the U - B and B - V val¬ 
ues suggest a KOIII classification rather than the G5 IV spectral 
type p reviously reported. Using spectroscopic data, iBopp et all 
(1 19931) confirmed the KOIII spectral type and determined an or¬ 
bital period of Rorb = 23.97 days and v sin / = 17 ± 2 km s“^. A 
record light-curve amplitude of 0 ?63 in V and 0^17 'mV - 1 
in January 1998 was reported by IStrassmei^ (Il999l) . At this 
state of high activity, spectroscopic observations were obtained 
at Kitt Peak National Obs ervatory (KPNO) that yielded the first 
Doppler image of XX Tri (IStrassmeiedI 19991) . It showed a gigan¬ 
tic, cool, high-latitude spot of elliptical shape with a temperature 
of around 1300 K below the photospheric temperature. Beside 
this superspot, a smaller, cool spot with a temperature of around 
4000 K and a warm equatorial spot on the adjacent hemisphere 
were reconstructe d. A 28-yr data set of photometric observations 
(lOlah et al.ll2014l). including about 20 years of APT photometry 
(see Istrassmeier et lH 1 1 9971 for details about APT), revealed a 
long-term modulation in V of 1TO5 from the deepest minimum 
to the overall maximum with a length comparable to the length 
of the data set. 

Most spotted stars investigated so far show solar-like differ- 
ential rotation, i.e., the equator rot ates faster than the pole (e.g., 
iBarnes et~^ 120051: iReinersI [20061) . In the last decade, several 
stars with antisolar differential rotation were detected (see, e.g., 
iKovari et~^ l2015l and references therein). As the differential 
rotation holds important information about the stellar dynamo 
working beneath the surface, theoretical developments need con¬ 
tinuous feedback from observations. It is particularly important 
to measure stars of different types because it is still not fully un¬ 
derstood how stellar dynamos work. In close binaries, such as the 
RS CVn-systems, tidal effects are thought to play an important 
role, as they help maintain the fast rotation and magnetic activity 
on high levels (cf.rScharlemannlll98lL[l982l:ISchriiver & Zwaanl 

ll99UlHolzwarth & Schussleill2002l). 


In this paper, we present a series of starspot distributions with 
different morphology, such as spot fragmentation and spot merg¬ 
ing, and with apparently a large range of variability timescales. 
We found evidence for starspot decay and active longitudes, and 
we estimate a magnetic cycle timescale. The paper is structured 
as follows. Sect, m describes the STELLA observations. Sect. [3] 
redetermines the absolute astrophysical parameters of XX Tri. 
Among these parameters is a revised orbit for the binary system, 
revised atmospheric abundances, and better constrained values 
for mass, radius, and age. Sect. |4] presents our Doppler imagery 
of a total of 36 maps over the course of six years. In this sec¬ 
tion, firstly we describe our Doppler-imaging code iMap and the 
parallel line-profile denoising algorithm. Secondly, we present 
and apply a phase-refilling scheme that uses data from the pre¬ 
vious or the following stellar rotation. Next, we present several 
tests of this scheme with partially artificial, partially real data. 
It includes inversion tests with phase gaps like in the real data, 
e.g., due to bad weather. Thirdly, we define the term spot area on 
the basis of photometric image analysis that fits a spot model to 
each Doppler image. These images are then analyzed in detail in 
Sect. [5] and discussed in Sect. (6] 

2. Spectroscopic observations 

Time-series, high-resolution echelle spectroscopy of XX Tri was 
taken with the L2-m STELLA telescopes and the STELLA 
Echelle Spectrograph (SES) on a nightly basis between July 
2006 and April 2012. The STELLA observatory is located at the 
Izana Observatory on Tenerife and operates fully ro botic with 
two L2-m telescopes (IStrassmeier et al.ll2004l2010bl) . The SES 
is a white-pupil spectrograph with an R2 grating with two off- 
axis collimators, a prism cross disperser and a folded Schmidt 
camera with an e2v 2kx2k CCD as the detector, the latter two 
items were replaced by a fully refractive camera and an e2v 
4kx4k CCD in mid-2012. In 2010, the SES fiber was moved to 
the prime focus of the second STELLA telescope (STELLA-II), 
while STELLA-I now hosts the Wide-Eield STELLA Imaging 
Photometer (WiESIP). E urther details o f the p erfor mance of the 
system were reported by ICranzer et al.l (l2010l) and IWeber et ^ 
(1201^ . 

We obtained a total of 667 usable spectra from six obser¬ 
vational seasons. Spectra cover the wavelength range from 388- 
882 nm with increasing inter-order gaps near the red end starting 
at 734 nm toward 882 nm before the camera and CCD exchange. 
The resolving power is R = 55,000 corresponding to a spectral 
resolution of 0.12 A at 650 nm. 

We set the integration time to 7200 s because of the relative 
faintness of the target for a Im-class telescope. Depending on 
weather conditions, the averaged signal-to-noise (S/N) ratios are 
between 50-300:1 for each spectrum, but typically 150:1. The 
SES spectra are automatically reduce d using the IRAE-b ased 
STELLA-SES data-reduction pipeline (IWeber et al.ll2008l) . We 
corrected the images for bad pixels and cosmic-ray impacts. 
We removed bias levels by subtracting the average overscan 
from each image followed by the subtraction of the mean of 
the (already overscan subtracted) master bias frame. We flat¬ 
tened the target spectra with a nightly master fiat, which itself 
is constructed from around 50 individual fiats observed dur¬ 
ing dust, dawn, and around midnight. After removal of scat¬ 
tered light, the one-dimensional spectra were extracted using an 
optimal-extraction algorithm. We then removed the blaze func¬ 
tion from the target spectra, followed by a wavelength calibra¬ 
tion using consecutively recorded Th-Ar spectra. Einally, the ex¬ 
tracted spectral orders were continuum normalized by dividing 


Article number, page 3of|25] 






























































A&A proofs: manuscript no. 25687_am_ok 


Table 1: Astrophysical properties of XX Tri (HD 12545). 


Parameter 

Value 

Based on 

Classification, MK 

KOIII 

Strassmeier (1999) 

Distance, pc 

160+32 

van Leeuwen (2001) 

Vmax, mag 

7.76 

Olah et al. 12014) 

V - Ic, mag 

1.18 

Olah et al. 12014) 

Rotation period, d 

;^24.0 

Strassmeier (1999) 

Orbital period, d 

23.96741 

radial velocities 

Inclination, deg 

60+10 

Strassmeier (1999) 

V sin /, km s“^ 

19.9 + 0.7 

spectrum synthesis 

Temperature, K 

4,620 + 30 

spectrum synthesis 

Log gravity, cgs 

2.82 + 0.04 

spectrum synthesis 

Metallicity, [Fe/H]© 

-0.13 + 0.04 

spectrum synthesis 

Micro turb., km s“^ 

1.5 

spectrum synthesis 

Macro turb., km s“^ 

3.0 

spectrum synthesis 

Radius, R© 

10.9+1.2 

from R sin i and i 

Luminosity, L© 

30!13 

from Mboi 

Mass, M© 

1.26 + 0.15 

evolutionary tracks 

Age, Gyr 

7.7 + 3.1 

evolutionary tracks 




-0.08 
^ - 0.10 
£ - 0.12 
-0.14 

I -0.16 

E -0.18 
- 0.20 

4000 4500 5000 5500 6000 4000 4500 5000 5500 6000 

HJD 2450000+ HJD 2450000+ 




Fig. 2: Stellar parameter determinations with PARSES. Shown 
are the obtained values for effective temperature, log gravity, 
metallicity, and v sin i for all spectra excluding 3-cr outliers. 


Notes. Values not cited in the third column were obtained in this pa¬ 
per. Note that the errors for spectrum-synthesis related parameters are 
internal errors. ^ + 0.0005, see text. 

with a flux-normalized synthetic spectrum of the same spectral 
classification as XX Tri. 

Fig.ia shows an example spectrum (Echelle order #89) from 
HJD 2,454,146.37 while Fig. [TJ) gives an overview of the time 
and phase sampling of all 667 STELLA spectra. 


3. Astrophysical parameters of XX Tri 

The revised reduction of the Hipparcos data (Ivan LeeuwenI 
l2007l) yielded a parallax of 6.24+ 1.02 mas and fixed the dis¬ 
tance of XX Tri (HIP 9630) to 160 pc. With an a pparent max¬ 
imum visual magnitude of 7^76 (lOlah et al.ll2014l) . the absolute 
visual magnitude of XX Tri is My = 1 ^ 58 ^q 4 q. W e took inter- 
stelkr absorption i nto account with () ^1 per 100 pc (IStrassmeied 
Il999l) . We note that lOlah et al.l (l2014l) estimated the reddening of 
XX Tri from all-sky infrared imaging and found a color excess 
E{B - V) of ^ 0^05, which leads to the same value of extinc¬ 
tion a s above. With a bolometric correction of -0^517 (iFlowed 
Il996l) . the bolometric magnitude of XX Tri is 1?06 and, with an 
absolute magnitude for the Sun of Mboi,© = 4^75, the luminosity 
must be approximately 30^g^ L©. 

Atmospheric surface stellar parameters (effective tempera¬ 
ture, log gravity, metalli city, and vsin/) were determined with 
the p rogram PARSES (lAllende Prietol |2004 iJovanovic et al.l 
l2013h . which is included in the SES data reduction pipeline. It 
fits synthetic spectra to a defined spectral region, in our case most 
echelle o rders between 480-750 nm, using MARCS model atmo¬ 
spheres (iGustafsson et al.l[20081) . We verified this approach b y 
applying it to the ELODIE library (iPrugniel & Soubir^l200ll) . 
and used linear regressions to the offsets with respect to the liter¬ 
ature values to correct the zero point of our PARSES results. In 
Fig. [2] these corrected values are shown for all spectra (exclud¬ 
ing 3-(T outliers). For more d etails of this procedure, we refe r 
to previous applications (e.g., IStrassmeier et all [2010aL l2012h . 
The mean values are Teff of 4620 ±30 K, a gravity log g of 
2.82 + 0.04, a vsini of 19.9 + 0.7 kms“\ and a metallicity of 
-0.13 + 0.04 dex relative to the Sun. Its errors are internal errors 


based on the rms of the entir e time series. The above effective 
temperature is close to what lOlah et aP (l2014l) obtained from 
VI photometry during maximum photometric brightness, but the 
metallicity differs by a factor of two. 

Fig. [2| also shows that the plotted values are systematically 
variable. In particular the effective temperature and line broad¬ 
ening clearly vary with the rotational period of the star and from 
season to season. There may be a similar trend in gravity, i.e., an 
increase during the first three observing seasons and a decrease 
during the last three seasons. A possible explanation could be re¬ 
lated to the assumption that a higher surface temperature means 
fewer cool spots, which implies a weaker internal magnetic field 
(and therefore magnetic pressure), and hence that the star con¬ 
tracts a bit and thus the surface gravity increases. We will in¬ 
vestigate this behavior in a forthcoming paper in more detail 
where we compare the observed broadband light curves with the 
photometric predictions from our Doppler maps along with the 
spectrum-integrated values from PARSES. 

To determine the mass and age of XX Tri, a trilinear inter- 
polation between stellar evolutionary t racks (iBertelli et al.ll2008l) 
based on a Monte Carlo (MC) method (lKunstled2008l) was used. 
Within the three-dimensional space (L, Teff, [Fe/H]) 10,000 ran¬ 
dom positions were generated, taking Gaussian errors into ac¬ 
count. For each generated position, we calculated the mass and 
age. The obtained mean values are a mass of 1.26 + 0.15 M© and 
an age of 7.7 + 3.1 Gyrs. Fig. [3] shows the evolutionary tracks 
interpolated to the metallicity of XX Tri including the star’s po¬ 
sition. 

A rotation period of 24.0 days together with a pro¬ 
jected rotational velocity of vsini = 19.9 + 0.7 kms“^ yields 
a minimum radius o f R sin i = 9.4 + 0.3 R©. With an inclina¬ 
tion of / 60+ 10° (IStrassmeied 1 19991) . the stellar radius is 

R = 10.9+1.2 R©. The unprojected equatorial rotational velocity 
would then be Veq = 23.0 km s“^ 

We use preliminary revised orbital elements from our 
radial-velocity fit of the STELLA data, see Fig. [T]:: 
Porb = 23.9674 + 0.0005 days, y = "25.389 + 0.031 kms-\ 
16.772 + 0.044 kms-\ a sin/ = 5.528+ 0.014x 10^ km, 
e = 0 (adopted), and f(M) = 0.0117 + 0.0001. Phase is always 
computed from a time of maximum positive radial velocity with 
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log T,„ 

Fig. 3: S telly evo lutionary tracks of low-mass stars from 
iBertelli et~^ (l2008l) together with the position of XX Tri. Each 
track is interpolated to a metallicity of -0.13 and covers the 
range from the ZAMS up to the RGB. The stellar mass (in M©) 
is indicated at the initial point of evolution. 

Table 2: Spectral lines used in the inversion process. 


Ion 

^(A) 

Ion 

^(A) 

Fei 

5049.820 

Fei 

5576.089 

Cut 

5105.537 

Cai 

5581.965 

Fei 

5198.711 

Cai 

5601.277 

Fei 

5232.940 

Fei 

6020.169 

Fei 

5302.300 

Fei 

6024.058 

Fei 

5307.361 

Fei 

6065.482 

Fei 

5324.179 

Cai 

6122.217 

Cri 

5345.796 

Fei 

6173.334 

Cri 

5348.315 

Fei 

6219.281 

Fei 

5367.466 

Fei 

6254.258 

Fei 

5383.369 

Fei 

6265.132 

Fei 

5393.167 

Fei 

6322.685 

Mni 

5394.677 

Fei 

6393.600 

Mni 

5420.355 

Fei 

6408.018 

Fei 

5434.524 

Fei 

6411.648 

Fei 

5445.042 

Fei 

6421.350 

Fei 

5497.516 

Fei 

6430.845 

Fei 

5501.465 

Cai 

6439.075 

Fei 

5506.779 

Cai 

6717.681 

Fei 

5569.618 

Fei 

6750.152 


the revised orbital period, 

HJD = 2,453,926.6663 + Ex 23.9674 . (1) 

Final orbital elements will be presented in our forthcoming paper 
where we remove the radial-velocity jitter due to spots and add 
three more years of data. 

The revised mass function, together with the primary mass 
of 1.26 M© and an orbital inclination of i ^ 60°, suggests a low- 
mass secondary star with a mass of 0.36 M©. Because the sec¬ 
ondary is not seen in the spectrum, the most likely secondary 
star is then a red dwarf of spectral type M. The most important 
astrophysical properties of XX Tri are summarized in Table [T] 


4. Doppler Imaging 

The spectral resolution of 55,000 combined with the vsini of 
just 20 km s“^ and the relative faintness of the star for a 1.2-m 
telescope, places XX Tri close to our limit for Doppler imaging. 
Significant effort is thus put to remove instrumental noise and 
other systematics from the data and to denoise the line profiles 
to prepare for the inversion. 


4.1. The iMap code 

All maps were computed with our Dopple r Imaging (Dl) and 
(Zeeman Doppler Imaging) ZDI-code iMap (ICarroll et al.ll2007l 
l2008l l2009l l2012h . Here w e give just a brief dp cription of the 
code, for further details see lCarroll et al.l (l2012h . The code per¬ 
forms a multiline inversion of a large number of photo spheric 
line profiles simultaneously. For the local line profile calcula- 
tion, the code utiliz es a full (polarized) radiative transfer solver 
(ICarroll et al.ll2008l) . The atomic line p arameters are taken from 
the VAFD dat abase (iKupka et H] 1999 ). We used Kurucz model 
atmospheres (ICastelli & Kuruczl 2004 ) which are interpolated 
for each desired temperature, gravity, and metallicity during the 
course of the inversion. Additional input parameters are v sin /, 
micro- and macroturbulence. 

Because of the typical ill-posed nature of the problem, 
an iterative regular ization based on a Landweber algorithm 
(ICarroll et al.ll2012h is implemented, having the advantage that 
no additional constraints are imposed in the image reconstruc¬ 
tion. For all temperature maps, the surface segmentation is set 
to a 5° X 5° equal-degree partition, resulting in 2592 segments. 
Because of the inclination of 60° a total of 432 segments are hid¬ 
den and therefore only 2160 segments are included during the 
inversion process. The code calculates the full radiative transfer 
of all involved line profiles for each surface segment depend¬ 
ing on the current effective temperature and atmospheric model. 
The surface temperature of each segment is adjusted according 
to the local (temperature) gradient information. The line profile 
discrepancy is reduced until a minimumis obtained. 


4.2. Line profile denoising 

We included 40 well-defined absorption lines simultaneously in 
our inversion, which are listed in Tabled These lines were cho¬ 
sen individually by investigating the stellar spectra and VAFD 
database and several other criteria, such as having a minimum 
line depth of 0.75 ///c, being almost blend-free, and having a 
good continuum stratification above 0.9 ///c- Additionally, all 
blends within + 1 A of each extracted line profile and a min¬ 
imum line depth of 0.1 are included in the inversion. As we 
have to deal with relatively low S/ N ratios, a wavele t analy¬ 
sis based on the a trous-Silgovithm (IStarck et al .11 19981 chapter 
1.4. 4 and references the rein) is implemented for further denois¬ 
ing. IStarck et"^ (1 19971) showed that for noisy data the wavelet 
transform is a powerful signal processing technique for spectral 
analysis. Each line profile, in our case the mean profile out of the 
40 individual lines, is split into so-called wavelet scales wj and 
a smoothed array Cp, whereas their sum represents the original 
spectrum co(A) = Cp(A) -h Xj=i For oach wavelet scale the 

standard deviation is determined and only signals above 3 cr are 
overtaken in the recomposition of the spectral line. 
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Fig. 4: Test #1 of the influence of phase gaps and filled gaps. 
Each image is shown in four spherical projections 90° apart, a) 
Input map. This map is identical to the first reconstruction from 
season 2007/08 (DI #8 from 2007.67) with a total of 21 phases, 
b) The reconstruction when ignoring phases no. 4-9 in the in¬ 
version process, c) The (absolute) difference a - b. d) The re¬ 
construction when the phase gap is filled with phases from the 
following stellar rotation, e) The (absolute) difference a- d. 


Fig. 5: Test #2 of the influence of phase gaps and filled gaps. 
Each image is shown in four spherical projections 90° apart, a) 
Input map. This map is identical to the first reconstruction from 
season 2007/08 (DI #8 from 2007.67) with a total of 21 phases, 
b) The reconstruction when ignoring phases no. 10-15 in the in¬ 
version process, c) The (absolute) difference a - b. d) The re¬ 
construction when the phase gap is filled with phases from the 
following stellar rotation, e) The (absolute) difference a- d. 


4.3. Phase selection and gap filling 

To later quantify the continuous evolution of spots based on con¬ 
secutive Doppler maps, we first deal with the inherent limita¬ 
tions of our phase coverage. In certain circumstances it is diffi¬ 
cult to compare consecutive maps that had different phase cov¬ 
erage and that even contained some larger observational gaps 
(several tenths of a phase) at different rotational phases at dif¬ 
ferent times. The effect of phase gaps on the recovery of in- 
divid ual spots had been simula ted by many authors in the past 
fe.g. jRice & Strassmeied[200Q[ and references therein). Gener¬ 
ally, Doppler imaging is very robust against small phase gaps 
but large phase gaps may introduce spurious spots at surface lo¬ 
cations not covered by the data. 

Fig -in and [5] show our simulations with iMap based on real 
data of XX Tri. During the season 2007/08, STELLA has cov¬ 
ered two consecutive stellar rotations completely with one obser¬ 
vation per night (DI #8 from 2007.67 and DI #9 from 2007.73, 


amounting to 21 phases from 23 nights per rotation). From 
DI #8, we removed six consecutive phases to create an artifi¬ 
cial phase gap of 90° (0^25), and compared the resulting map 
with the original map (Fig. |4^-c). The two darkest and biggest 
spots at phases around 17 p 25 and 17 p 55 could not be separated 
anymore. The larger spot loses a big part of its area, which is 
seen in the difference map in Fig. with a temperature simi¬ 
lar to the difference between photospheric and spot temperature. 
In the next step, we filled these gaps with observations from the 
following stellar rotation (DI #9) and again compared the result¬ 
ing map with the original map (Fig. |4^-e). All individual spots 
are now reconstructed with no changes of their size or temper¬ 
ature exceeding the expected errors driven by the S/N of the 
data. Fig. [S^-c shows another simulation of the same data with 
an artificial phase gap of 90° but at a different rotational phase. 
Here, the smaller spot at phase around 17 p 55 has almost com¬ 
pletely vanished, whereas it is recovered for the case with gap 
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filling (Fig.[5]i-e). We further verified this method on a few other 
Doppler images with the same result (not shown). 

We conclude that large 90°) phase gaps have a strong im¬ 
pact on the recovery of the global stellar spot distribution. It af¬ 
fects not only their size and shape but also their location. For 
smaller spots located within or near the missing phases it af¬ 
fects even their existence. However, when we compare the orig¬ 
inal map with the maps where the missing phases were filled by 
phases from one rotation earlier or later, we see a significantly 
better agreement. This is the case even if the spot distribution had 
evolved in the meantime. Based on these simulations, we have 
a simple method to minimize the systematics due to phase gaps. 
This is only a first-order approximation because the spot evolu¬ 
tion within the missing surface area remains unknown. Fig. 1C. II 
shows the phase coverage and gap filling of each Doppler image 
for all seasons, if applicable. Real phase gaps ranged between 
1-13 d or 0 p04-0p 54 at the extremes, but with typical values more 
near 4 d or OH?. 

For two out of the 36 Doppler images, we had to use phases 
from two rotations earlier and/or later, as there were no observa¬ 
tions closer in time available. In the first case (DI #5) we filled 
a gap of 5 d (0 p 21) with two phases. To minimize a possible 
smearing effect, one phase was taken from two rotations earlier 
and one from two rotations later. In the case of DI #20, we had 
to deal with a small number of existing phases in addition. In 
general, we would not have taken this rotation into account, but 
in this particular case, we would have only four maps for the ob¬ 
serving season 2007/08, two consecutive maps at the beginning 
and two consecutive maps at the end of the season, which would 
have severely limited our coverage for this season. We will redis¬ 
cuss this explicitly in Sect. l5.ll in terms of spot-area evolution. 

Table [3] summarizes our Doppler-image log. The year indi¬ 
cates the mid-time of the Doppler image, followed by its Helio¬ 
centric Julian Date (HJD) range, the time At elapsed in days, the 
number of spectra N, the largest phase gap in number of con¬ 
secutive days (the rotation period is 24 d), and the number of 
observations from following or preceding stellar rotations that 
were used to fill phase gaps. In total, we obtained 36 individual 
Doppler images, each with between 11 to 24 spectra. 

4 . 4 . Image analysis: definition of the spot area 

One parameter to be extracted from each image is the surface 
area of a spot. Its measurement depends on the definition where a 
spot ends and where the undisturbed photosphere begins. Our in¬ 
version code is completely free in the reconstruction and usually 
recovers an irregular spot morphology. Furthermore, the images 
contain small-scale structures, which seem to appear or disap¬ 
pear from one rotation to another and/or move to different longi¬ 
tudes/latitudes on timescales that are likely shorter than a stellar 
rotation. These structures often appear as hot spots or as a pair 
of a hot and a cool spot (but not necessarily along the same iso- 
radial line which would indicate an artifact). The appearance of 
elongated appendages to the polar spot, or spots that appear con¬ 
nected to another spot, are further complications for determin¬ 
ing a spot’s area. Because of inherent limitations of the spatial 
surface resolution due to the low v sin /, as well as to a lesser 
extent also due to the choice of regularization during the inver¬ 
sion, these spot configurations cannot easily be separated from 
each other anymore. Simple area integration of the disk within a 
certain temperature difference would thus lead to erroneous spot 
areas. 

We define the spot area by fitting artificial spot models to 
the maps. Our artificial spots have circular shapes of arbitrary 


size, but a constant temperature of 3500 K. From photometric 
spot-modeling with spot temperature as free parameter the de¬ 
rived temperature difference AT = Tphot - T'spot were determined 
to 11 on K d^rassmeier & Olahll 992h and 1 280 K dRkeH^l 995h as 
well as in the range of 650-1200 K (iHampton et al.l 19961) for var- 
ious epochs. Furthermore, the superspot on t he first Doppler ini - 
age had a temperature difference of 1300 K (IStrassmeiedI 19991) . 
These values are in agreement with the derived spot temperatures 
of ^ 3500 K from our Doppler images for long-lived spot struc¬ 
tures. As our focus lies in the evolution of starspots, i.e., their 
decay or growth, we investigate mainly large-scale spot struc¬ 
tures that are repeatedly reconstructed from one stellar rotation 
to the next. 

Our method is based on the spot-modeling procedure used 
in light-curve analysis fe.g.. iRibarik et al .1120031) . where an ap¬ 
propriate number of spots with circular shape and a defined tem¬ 
perature is taken as input. An initial guess of the spot’s location 
and radius was taken directly from the observed Doppler images. 
With these starting values, we calculated the best fit with an area- 
and temperature-weighted Monte Carlo (MC) method and thus 
extracted a definition-dependent, best-effort spot location and 
area. Within the three-parameter space (longitude, latitude, ra¬ 
dius) 10,000 random positions were generated, using a range of 
15° for each parameter, and then cross correlated with the orig¬ 
inal Doppler map. From the best 100 correlation maps (which 


corresponds to ^^N), the mean values and their standard devia¬ 
tions are determined. In the following, spot area always refers to 
the spot area deduced from this analysis. 

The area of the individual spots from the spot-model fits is 
summarized in Table [3] in units of solar hemispheres (SH). To 
estimate the quality of the spot-model fit, we compared the to¬ 
tal spotted area between the spot model and the Doppler image. 
To determine the total spotted area of the Doppler image, a tem¬ 
perature weighting by analogy to the MC method was used. We 
obtain the total spotted area such that each spotted segment i at¬ 
tributes a certain fraction of its area at to the total spotted area 
dependent on its temperature given by 



i 


Gi 


(TeS - Tj) 

(^^)max 


( 2 ) 


where (Ar)niax = (T^phot - 7^spot)max = 1120 K. The small-scale 
surface structures, cool and hot spots, were not counted into the 
total spotted area. The quality of the spot-model fits is given in 
Table [3] in terms of the uncertainty of the spot models, where all 
our spot models lie within 1 cr. The absolute scale of the area 
units in m^ is set by the stellar radius of 10.9 R©. Formally, the 
stellar surface of XX Tri is 724 Gm^ and thus 118.8 times the 
surface of the Sun or 1 SH 0.8 % of an hemisphere of XX Tri. 


4 . 5 . Doppler images and spot models 

Fig.[6]is a representative figure for our results and shows all maps 
for the observing season 2006/07. Seven consecutive Doppler 
images are reconstructed from a total of 9.5 stellar rota¬ 
tions. All maps show a large polar spot with a temperature of 
3500 K. During this season, the polar spot drifted apart and 
changed its morphology from almost circular to an elongated 
spot form. This drift could be a si gn of differential rotation and 
is investigated further in Section 15.31 Furthermore, a smaller 
high-latitude spot with a temperature of around 3800 K is re¬ 
constructed. It is seen that the larger spot approached the smaller 
spot. Because of the inherent technical limitations of Doppler 
imaging, the surface resolution near the rotational pole is poor 


Article number, page 7 of|25] 













A&A proofs: manuscript no. 25687_am_ok 


(a) 


(p= 0.34 <p= 0.59 


(C) 


(d) 


(p= 3.35 (p= 3.6 


(e) 


(f) 


ip= 6.89 <p= 7.14 


(g) 


0.84 (p= 1.09 




2.85 (p= 3.10 


3.85 (p= 4.10 


5.85 (p= 6.10 


7.39 (p= 7.64 


9.15 <!>= 9.40 


I 

I 


5000 „ 
4625 ^ 

3 

4250 « 
o 

3875 E 

CD 

3500 


5000 ^ 
4625 ‘o 

3 

4250 « 
a> 

3875 E 

CD 

3500 










3875 E 


Fig. 6: Doppler images of XX Tri for the observing season 
2006/07. Each image is shown in four spherical projections sep¬ 
arated by 90°. The rotational shift between consecutive images is 
corrected, i.e., the stellar orientation remains the same from map 
to map and from season to season. The time difference between 
each Doppler image is indicated in units of rotational phase 0. 
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Fig. 7: Spot-model fits of the Doppler images in Fig.[6l Each spot 
is shown with different color/contrast for better visualization. 


and thus spot separations not well constrained. Therefore, one 
cannot say whether the large spot is a monolithic structure or 
being a conglomerate of several smaller spots. Besides the po¬ 
lar spot, scattered small cool and/or hot spots are visible at lat¬ 
itudes between 0-60° with absolute temperatures between 4200 
to 5000 K. 
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Table 3: Doppler image log, fit quality, surface temperature, and spot areas. Detailed information is given in Section l43ll4.5l 


DI 

# 

Year 

HJD range 
(2450000+) 

Ar 

(d) 

N 

0-gap 

(d) 

N-mr 

RMS'’ 

RMS" 

(t^spot) 

Temperature (K) 

Tmax Tniean 

Spot A 

Spot B 

Spot area (SH)"^ 
Spot C Spot D 

Spot E 

Spot Total 

1 

2006.58 

3935-3963 

28 

15 

4 

2 

0.0048 

0.6 

4667 

4470 

9.2 ±0.3 

1.1 +0.6 

1.3+ 0.8 

— 

— 

11.5 + 1.0 

2 

2006.64 

3959-3980 

21 

19 

3 

0 

0.0046 

0.6 

4731 

4477 

8.5+ 0.5 

1.5+ 0.5 

1.2+ 0.7 

— 

— 

11.2 + 0.9 

3 

2006.71 

3966-4015 

49 

14 

5 

3 

0.0047 

0.7 

4796 

4489 

8.1 + 0.6 

2.1 +0.6 

1.6+ 0.9 

— 

— 

11.8 + 1.0 

4 

2006.78 

4007-4050 

43 

18 

7 

3 

0.0047 

0.7 

4961 

4476 

7.2+ 0.7 

2.6+ 0.5 

2.4 + 0.7 

— 

— 

12.2 + 0.8 

5 

2006.91 

4011-4108 

9r 

17 

5 

3 

0.0047 

0.8 

5059 

4477 

6.0+ 0.5 

3.3+ 0.5 

3.0+ 1.1 

— 

— 

12.3 + 1.1 

6 

2007.01 

4092-4138 

46 

16 

5 

2 

0.0044 

0.5 

5052 

4492 

5.4+ 0.6 

3.7+ 0.6 

3.1 + 1.2 

— 

— 

12.2+1.3 

7 

2007.14 

4124-4161 

37 

19 

5 

2 

0.0043 

0.7 

5125 

4486 

5.4+ 0.7 

3.9+ 0.6 

3.2+ 1.1 

— 

— 

12.5 + 1.3 

8 

2007.67 

4333-4356 

23 

21 

3 

0 

0.0040 

0.4 

5129 

4541 

5.5+ 0.5 

1.6+ 0.4 

— 

— 

— 

7.1 +0.8 

9 

2007.73 

4357-4380 

23 

21 

2 

0 

0.0035 

0.4 

5027 

4537 

5.2+ 0.4 

1.8+ 0.7 

— 

— 

— 

7.0+ 0.9 

10 

2007.87 

4362-4475 

113" 

15 

13 

8 

0.0033 

0.7 

4912 

4538 

4.1 + 0.7 

2.6+ 0.3 

— 

— 

— 

6.7 + 0.7 

11 

2008.05 

4468-4512 

44 

19 

4 

1 

0.0037 

0.6 

4685 

4534 

4.3 + 0.7 

2.5 + 0.4 

— 

— 

— 

6.8+ 0.7 

12 

2008.12 

4481-4516 

35 

15 

5 

2 

0.0038 

0.6 

5188 

4524 

4.1 + 0.5 

2.6+ 0.5 

— 

— 

— 

6.6 + 0.4 

13 

2008.53 

4649-4671 

22 

20 

2 

0 

0.0040 

0.5 

4903 

4503 

4.8+ 0.2 

3.2+ 0.3 

— 

— 

— 

8.0 + 0.4 

14 

2008.60 

4673-4696 

23 

22 

3 

0 

0.0037 

0.5 

4839 

4501 

5.5+ 0.3 

3.0+ 0.4 

— 

— 

— 

8.5+ 0.5 

15 

2008.73 

4722-4744 

22 

14 

5 

0 

0.0038 

0.7 

4845 

4501 

6.4 + 0.4 

2.0+ 0.2 

— 

— 

— 

8.4 + 0.4 

16 

2008.80 

4726-4790 

64 

15 

7 

5 

0.0034 

0.5 

4745 

4497 

7.5+ 0.5 

1.6+ 0.4 

— 

— 

— 

9.1 +0.6 

17 

2008.88 

4775-4798 

23 

16 

3 

0 

0.0035 

0.4 

5057 

4503 

8.0+ 0.7 

0.7 + 0.2 

— 

— 

— 

8.7+ 0.8 

18 

2008.94 

4799-4822 

23 

15 

4 

0 

0.0032 

0.3 

5004 

4508 

8.1 + 0.6 

0.6+ 0.3 

— 

— 

— 

8.8+ 0.7 

19 

2009.07 

4844-4873 

29 

16 

5 

2 

0.0035 

0.9 

5193 

4493 

8.2+ 0.4 

1.2+ 0.3 

— 

— 

— 

9.4 + 0.5 

20 

2009.60 

5039-5066 

27 

18 

5 

1 

0.0033 

0.3 

4918 

4522 

7.8+ 1.2 

2.9 + 0.4 

— 

— 

— 

10.7 + 1.3 

21 

2009.67 

5039-5083 

44 

17 

6 

2 

0.0035 

0.7 

4825 

4500 

8.1 + 1.0 

3.3+ 0.4 

— 

— 

— 

11.3 + 1.1 

22 

2009.86 

5135-5161 

26 

17 

4 

1 

0.0039 

0.9 

5061 

4495 

8.3 + 1.1 

4.7 + 0.5 

— 

— 

— 

13.0+1.1 

23 

2009.93 

5154-5173 

19 

16 

11 

3 

0.0036 

0.4 

5001 

4511 

7.4+ 1.1 

4.8+ 0.3 

— 

— 

— 

12.2+1.1 

24 

2010.06 

5207-5224 

17 

17 

7 

0 

0.0035 

0.6 

5146 

4506 

6.8+ 0.7 

4.9 + 0.6 

— 

— 

— 

11.6+1.0 

25 

2010.59 

5401-5447 

46 

12 

6 

4 

0.0029 

0.6 

4935 

4535 

4.5 + 0.7 

3.3+ 0.3 

0.9 ± 0.6 

— 

— 

8.7+ 0.9 

26 

2010.67 

5429-5452 

23 

19 

2 

0 

0.0031 

0.4 

5088 

4538 

4.0+ 0.7 

3.1 +0.4 

1.1 ±0.6 

— 

— 

8.1 +0.9 

27 

2010.74 

5439-5496 

57 

13 

7 

3 

0.0029 

0.9 

4967 

4532 

3.2+ 0.4 

2.4 + 0.6 

1.5 ±0.7 

— 

— 

7.2+ 0.9 

28 

2010.84 

5467-5512 

45 

17 

7 

2 

0.0030 

0.9 

4943 

4522 

2.6 + 0.4 

1.5+ 0.6 

2.4 ± 0.6 

— 

— 

6.5 + 0.7 

29 

2011.12 

5595-5614 

19 

11 

5 

0 

0.0021 

0.5 

4620 

4513 

— 

— 

— 

5.1 +0.5 

2.3 + 0.4 

7.4+ 0.5 

30 

2011.61 

5774-5797 

23 

18 

3 

0 

0.0033 

0.3 

4758 

4493 

4.7 + 0.6 

2.9 + 0.4 

1.2+ 0.6 

— 

— 

8.8+ 0.7 

31 

2011.68 

5799-5822 

23 

24 

1 

0 

0.0032 

0.4 

4756 

4501 

4.5 + 0.7 

2.5+ 0.8 

1.7+ 0.6 

— 

— 

8.8+ 0.8 

32 

2011.81 

5830-5879 

49 

17 

7 

2 

0.0029 

0.8 

4649 

4509 

3.7+ 0.7 

1.4+ 0.7 

3.1 +0.4 

— 

— 

8.2+ 1.0 

33 

2011.88 

5871-5894 

23 

16 

4 

0 

0.0025 

0.9 

4722 

4501 

3.0+ 0.5 

2.7+ 0.8 

3.4+ 0.3 

— 

— 

9.1 +0.8 

34 

2012.01 

5919-5940 

21 

20 

3 

0 

0.0027 

0.4 

4664 

4486 

3.6+ 0.7 

3.4+ 0.9 

2.6 + 0.4 

— 

— 

9.6+ 1.0 

35 

2012.08 

5944-5967 

23 

18 

3 

0 

0.0026 

0.7 

4620 

4506 

3.3+ 0.7 

4.2 + 0.9 

1.7+ 0.4 

— 

— 

9.2+ 1.1 

36 

2012.19 

5983-6006 

23 

17 

3 

0 

0.0031 

0.3 

4767 

4477 

3.6+ 0.7 

5.8+ 0.6 

0.6 + 0.4 

— 

— 

10.0 + 0.8 


Notes. Number of phases/observations borrowed from the following and/or preceding stellar rotation. Deviation between observed and calculated line profiles. Deviation of spot area between 
spot models and Doppler image normalized by the standard deviation of the total spotted area. In units of solar hemispheres (SH). Including observations from two stellar rotations earlier 
and/or later. 
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Compari ng our new Dopple r images with the first image 
published by IStrassmei^ (Il999l) shows a very good agreement 
of spot locations and spot temperatures, despite 10-year time 
difference. In both cases, the largest spot appears near the pole 
with a temperature of around 3500 K. Smaller cool/hot spots 
appear on lower latitudes, again with comparable temperature 
differences. The super spot from 1998 with an area of approx¬ 
imately 11 % of the entire stellar surface is much larger than 
the largest recovered spots that we present. This is expected be¬ 
cause the star appeared to be in its brightest stage ever in 2009. 
In parallel, the rotationally-modulated photometric light curves 
show a small amplitude compared to that in 1998. During the 
season 1998/99 a record amplitude of 0% in V was observed 
in comparison to an amp litude of < 0?4 d uring the seasons from 
2006/07 to 2011/12 tsee lOlah et alllml . 

In Fig. |7] the best spot-model fits for observing season 
2006/07 are shown. Three artificial spots were sufficient to reach 
a high correlation for all Doppler images in this season. Two 
spots were utilized to represent the large polar spot and to model 
its drift during almost ten stellar rotations. The third spot repre¬ 
sented the smaller high-latitude spot on the opposite hemisphere. 
All other cool and/or hot spots at lower latitudes with no or only 
very short continuous appearance were ignored and therefore not 
implemented in the spot-model analysis. 

Doppler images and spot models for the observing seasons 
2007-12 are shown in Fig. lA.lllA.lOl The observed and inverted 
line profiles of each Doppler image are given in Fig. lB.ll and lB.2l 
Table[3]lists the fit quality between observed and calculated 
(xi^c) profiles. 


RMS = 




1 ” 

- 'y,(xi,o - , 

n ^ 


(3) 


where n is the total number of points of all line profiles used 
during the inversion. Table [3] also tabulates the maximum and 
mean surface temperatures of each Doppler image. 


5. Results 

5.1. Spot area evolution 

In each observing season, we see both decay and growth of 
individual spots. From these spots we infer a linear area de¬ 
cay/formation law of form 

dA(t)/dt = D. (4) 

As k nown from sunspot-decay studies (e.g., iMartmez Pilletl 
l2002h , the application of a linear law is the most appropriate way 
to describe sunspot decay as well as sunspot growth rates. Fig.[H] 
shows the area evolution of the total spotted area of XX Tri as 
well as of the individual spots from the spot models for each sea¬ 
son. The numerical values of D for each spot together with the 
overall mean are summarized in Table (U 

Season 2006107. Three artificial spots are sufficient to match 
the observed spot distribution in this season. The three spots 
appear partly merged and make up for the elongated polar- 
spot appendage in the Doppler images. In the case of overlap¬ 
ping spots, the overlapped area refers to the larger spot. The 
two overlapping spots, A and B, fragment, while the larger 
spot A shrinks and the smaller spot B waxes. Between DI #1- 
6 (At ^ 6.5 Frot) spot A loses around 40 % of its area from 
9.2 to 5.4 SH (D = -0.025 + 0.003 SH/day), while the smaller 
spot B increases to more than three times its area from 1.1 


to 3.7 SH (D = +0.017 + 0.004 SH/day). It seems that spot A 
“feeds” spot B, suggesting flux transport between these two 
spots. Between DI #6-7 both spots A and B remain almost 
constant in area. Spot C is located on the opposite hemisphere 
with a longitudinal shift of around 180°. Between DI #2-5 
(At ^5.5 Frot) it waxes to almost three times its area from 
1.2 to 3.0 SH (D = +0.020 + 0.012 SH/day). Afterwards, spot 
C remains almost constant in area. If we were to exclude 
DI #5 (see Sect, the determined values of D would be 
-0.024 + 0.004 SH/day for spot A, +0.017 + 0.005 SH/day for 
spot B, and +0.016 + 0.010 SH/day for spot C, respectively. 

Season 2007108. We concentrate on the two larger spots 
(at phases around 17 p 25 and 17 p 55) for a meaningful anal¬ 
ysis of spot area evolution. Both spots A and B are sepa¬ 
rated by around a quarter of rotation in DI #8. During the 
observing season they merge, as the larger spot A rotates 
much slower near the pole than the smaller spot B, which 
is located at mid-latitudes. This indicates that differential ro¬ 
tation is detectable as shown in Sect. lO Between DI #8- 
10 (At = 3 Frot) spot A loses around 25 % of its area from 
5.5 to 4.1 SH (D =-0.021 + 0.012 SH/day), while spot B 
increases to almost two times its area from 1.6 to 2.6 SH 
(D = +0.014 + 0.007 SH/day). This phenomenon of two spots 
interacting with each other, while one spot is decaying and the 
other is growing, has also been detecte d in t he previous season. 
If we were to exclude DI #10 (see Sect. 14.31) . a reliable determi¬ 
nation of D for this season would not have been possible because 
of the time sampling of the Doppler images. However, removing 
it from the entire time series has no impact on the mean decay 
rate. 

Season 2008109. Two spots are sufficient to characterize 
the spot evolution during this season. The two spots, A and 
B, are very close and appear connected to each other as seen 
in DI #13. The large spot A increases from 4.8 to 8.0 SH 
(D = +0.026 + 0.004 SH/day) between DI #13-17 (At ^ 5 F^ot) 
and afterwards remains almost constant in area. Within the same 
time span, the smaller spot B loses in area from 3.2 to 0.7 SH 
(D = -0.020 + 0.003 SH/day). After an almost complete decay, 
it starts to increase to two times its area from 0.6 to 1.2 SH 
(D = +0.013 + 0.009 SH/day) between DI #18-19 (At ^ 2 F^ot). 
As in the previous season, indications of differential rotation are 
seen. Spot A is located nearer to the pole than spot B and there¬ 
fore both spots get separated from each other. 

Season 2009110. Again, two spots are adequate to match the 
observed spot distribution throughout this season. Both spots 
A and B are located at opposite hemispheres with a longitu¬ 
dinal shift of around 180°. Between DI #20-22 (At = 4 Fj-ot) 
the smaller spot B waxes to around 160% its area from 2.9 
to 4.7 SH (D =+0.019 + 0.007 SH/day), whereas the larger 
spot A remains almost constant in area. Between DI #22-24 
(At = 3 Frot) the larger spot A loses in area from 8.3 to 6.8 SH 
(D = -0.019 + 0.016 SH/day), whereas the smaller spot B re¬ 
mains almost constant in area. 

Season 2010/11. During this season three spots are required 
to characterize the spot evolution. The two large polar spots, 
A and B, are located at opposite hemispheres with a longitu¬ 
dinal shift of around 150° and are moving toward each other. 
Between DI #25-28 (At ^ 4 Fj-ot) both spots A and B lose in 
area from 4.5 to 2.6 SH (D = -0.021 + 0.008 SH/day) and 3.3 
to 1.5 SH (D = -0.020 + 0.007 SH/day), respectively. During 
the same time span, the smaller spot C increases from 0.9 to 
2.4 SH (D = +0.017 + 0.009 SH/day) and moves toward higher 
latitudes. There is a large time gap of around four rotations be¬ 
tween DI #28-29. During this time span, the spot configura- 
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Fig. 8: Spot area evolution on XX Tri from 2006 to 2012. Shown are the seasonal evolution of the individual spots (dotted colored 
lines) for each observing season. The solid (colored) lines represent linear fits to the decay or growth of a spot. The black dotted 
line represents the total spotted area. The spot area is given in solar hemispheres on the left axis (1 SH = 3.05 Gm^) and relative to 
the total area of a stellar hemisphere of XX Tri on the right axis. 


tion obviously changes to a two-spot-model. As it is not clear 
whether one can identify these spots with spots from DI #25-28, 
we decide to regard them separately. These two polar spots, D 
and E, are located at opposite hemispheres with a longitudinal 
shift of around 180°. This spot configuration is almost identical 
with the first DI (#30) of the following observational season. 

Season 2011112. Again, three spots are sufficient to match the 
observed large-scale spot distribution during this season. As in 
the first season, the three spots appear partly merged and make 
up for the elongated polar-spot appendage in the Doppler im¬ 
ages. The two spots, A and B, are located at opposite hemi¬ 
spheres with a longitudinal shift of around 180°. This spot dis¬ 
tribution is very similar to that at the end of the previous season. 
The larger spot A fragments into two smaller spots, A and C, 
between DI #30-32. Spot C is located very close to the pole and 
rotates much slower than the other two spots A and B. Between 
the time of DI #33-36 spot C merges with spot B. Therefore, 
it again suggests flux transport between spot A and B. Spot A 
loses in area from 4.5 to 3.0 SH (D = -0.021 + 0.012 SH/day) 
between DI #31-33 (At ^ 3 Frot) and remains almost constant 
in size afterwards. Spot B loses in area from 2.9 to 1.4 SH 
(D = -0.021 ± 0.011 SH/day) between DI #30-32 (At = 3 Frot) 
and waxes between DI #32-36 (At ^ 6 Frot) up to 5.8 SH 
(D = +0.030 + 0.006 SH/day). Spot C increases in area up to 
3.4 SH (D = +0.023 + 0.006 SH/day) between DI #30-33 and 
decays between DI #33-36 (At ^ 4.5 Frot) almost completely 
(D = -0.024 + 0.005 SH/day). 

Finally, if we exclude Doppler images #5 and #10, our 
spot area evolution analysis would lead to an identical mean 
value for spot decay (<D> = -0.022 + 0.002 SH/day) and 


only a marginally increased mean value for spot formation 
(<D> = +0.022 + 0.002 SH/day). Therefore, we include them in 
our analysis. Furthermore, we repeated the entire analysis also 
for the scattered small cool and/or hot spots at low latitudes. 
Their respective values of D (for small cool spots +0.02 for 
growth and decay, respectively; hot spots +0.03) scatter within 
the range of the large-scale spots, but the time sampling is such 
that we can not determine a true beginning nor ending of the 
evolution. 

5.2. Active longitudes 

Active longitudes are longitudes on which spots occur pref¬ 
erentially. The analysis of long-term photometric observations 
as well as time-series Doppler imaging revealed active longi¬ 
tudes on several stars. iBerdvugina & T^uominenl (1 19981) found 
permanent active longitudes on four RS CVn stars. If two ac¬ 
tive longitudes, which are typically separated by 180°, change 
their spot activity a so-called “fiip-fiop” occured. This kind of 
a phenomenon w as first noticed on the late-type giant FK Com 
(jjetsu et al.lll99ll) . The average time between these phenomena 
is referred to as fiip-fiop cycle and has been observed to be in the 
range of a few years up to a d eca de. An observational overview 
is given in iBerdvuginal (l2007l) and iKorhonen & JarvinenI (l2007l) . 
In binary sta rs, a longitudinal dependence due to tidal effects 
is suggested (iHolzwarth & Schussledl2002b . Observations show 
that in binaries preferred longitude s exist on giant c omponents 
mostly at the substellar points fe.g. JOlah et al.11200^ . 

Fig.EIt-f shows the mean longitudinal distributions of all 
individual spots from our spot-models for each season. There is 
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Fig. 9: Active longitudes on XX Tri from 2006 to 2012. a) The overall mean distribution of the spot area (histogram in bins of 15°) 
is shown. The spot area is given in solar hemispheres on the left axis (1 SH = 3.05 Gm^) and relative to the total area of a stellar 
hemisphere of XX Tri on the right axis. The dashed line represents the phase toward the companion star, whereas the dotted line 
represents the phase in the opposite direction, b) Spot longitudes as a function of time. Spots A-E from all Doppler images are 
shown. The filled dots represent the larger spot at any given time. The dashed and dotted lines represent again the phases toward the 
companion star and in the opposite direction, respectively. 


Table 4: Individual spot decay and growth rates. Derived val¬ 
ues represent a linear area decay law of form dA{t)ldt = D. D is 
given in solar hemispheres (SH) per day (1 SH = 3.05 Gm^). 


Season 

Spot 

DI# 

D (SH/day) 

type 

2006/07 

A (blue) 

1-6 

-0.025 + 0.003 

decay 


B (red) 

1-6 

-rO.017 + 0.004 

growth 


C (green) 

2-5 

-rO.020 + 0.012 

growth 

2007/08 

A (blue) 

8-10 

-0.021+0.012 

decay 


B (green) 

8-10 

+0.014 + 0.007 

growth 

2008/09 

A (blue) 

13-17 

+0.026 + 0.004 

growth 


B (green) 

13-17 

-0.020 + 0.003 

decay 


B (green) 

18-19 

+0.013 + 0.009 

growth 

2009/10 

A (blue) 

22-24 

-0.019 + 0.016 

decay 


B (green) 

20-22 

+0.019 + 0.007 

growth 

2010/11 

A (blue) 

25-28 

-0.021+0.008 

decay 


B (green) 

25-28 

-0.020 + 0.007 

decay 


C (red) 

25-28 

+0.017 + 0.009 

growth 

2011/12 

A (blue) 

31-33 

-0.021+0.012 

decay 


B (green) 

30-32 

-0.021+0.011 

decay 


B (green) 

32-36 

+0.030 + 0.006 

growth 


C (red) 

30-33 

+0.023 + 0.006 

growth 


C (red) 

33-36 

-0.024 + 0.005 

decay 

2006-12 



<-0.022 + 0.002> 

decay 




<+0.021 + 0.002> 

growth 


clear evidence for preferred longitudes during each season but 
significantly spread out in location because of individual spot 
evolution. A particularly well-defined pair of active longitudes 
separated by 180° is seen in season 2009/10. Averaging the lon¬ 
gitudinal spot distribution from all 36 Doppler images, we find 
the most spotted longitude to appear on average in phase toward 
the unseen companion star (+ 90°; Fig.|9t). Fig.l^) plots the spot 
centers from our spot-model fits (spots A-E) as a function of 
time. The larger spot always appears in alternating hemispheres 


at locations either close to the phase toward the companion star 
or shifted by around 180°. The larger spot appears in between 
these phases near quadrature possibly only between DI #29-32. 
Nevertheless, we interpret this behavior as a flip-flop and esti¬ 
mate a tentative period of around two years. 

5.3. Differential rotation 

Tracking sunspots is a classic technique to measure solar differ¬ 
ential rotat ion and other surface velocity fields like meridional 
flows fe.g.. lWoh]|l2002t iBraisa et al.ll20()2h . In case of XX Tri a 
large number of temperature surface maps with unprecedented 
good sampling is available, and therefore may enable us to re¬ 
veal a similarly accurate differential rotation law by tracking in¬ 
dividual starspots. Differential rotation has been detected on a 
number of stars by cross-corre lating consecutive Doppler images 
in longitudinal direc tion fe.g.. lDonati & Collier Cameronll 19971: 
iKovari et~n]l2007l3) . Applying this method to our data, we re¬ 
construct between three and six cross-correlation-function (ccf) 
maps per observing season, which we average to increase their 
validity. We did not use the ccf map with DI #29, the last map 
in season 2010/11, because the time span to the previous map 
(#28) is comparable large, approximately four rotational periods. 
Within such a time span, we expect significant local spot evolu¬ 
tion. Fig. lD.ll shows the average ccf maps for each observational 
season. The resulting grand average ccf map, which consists of 
29 ccf maps in total, is given in Fig.fTOb. 

We determined the correlation peak for each longitudinal 
stripe of 5° width with a Gaussian profile and fitted a standard 
differential rotation law of the form 

Q(Z7) = Qeq-AQsin4fe), (5) 

which is usually used for differential rotation measurements on 
stars. The parameter ^(b) represents the angular velocity at lat¬ 
itude b, while AQ = Deq “ ^poie represents the difference be¬ 
tween the angular velocities at the equator and at the pole, re¬ 
spectively. The surface shear parameter a is defined as AQ/Deq, 
and the lap time as the reciprocal of the rotational shear, i.e., the 
time it takes for the equator to do a full lap more than the pole. 
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(a) 


5.4. Stellar cycle prediction 


Season 2006-2012 



Longitude shift (deg) 


(b) 


Season 2006-2012 
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Fig. 10: a) Grand average cross-correlation-function map from 
2006 to 2012. The map represents the average ccf map of all 
seasonal ccf maps from Fig. ID. 11 thus consisting of 29 ccf maps 
in total, b) Global differential rotation signature. The dots are 
the correlation peaks per 5°-latitude bin and their error bars are 
defined as the FWHM of the corresponding Gaussian fits. The 
dashed line represents a fit using Eq. [5l whereas the solid line 
represents a fit using Eq. [6j The parameters for each fit are sum¬ 
marized in Table [21 


Alternatively, we fit a differential rotation law of the form 


fl(b) = Qeq + ^1 -L Q 2 sm^(b ), 


( 6 ) 


which is usually used for differential rotation measurements on 
the Sun. In this case the angular velocity at the pole is defined as 
Qpoie = Deq -L Qi -L Q 2 - Fig. ID. 21 and Fig.fT^ show the observed 
differential rotation pattern determined from the ccf maps to¬ 
gether with the best fit of the differential rotation following Eq.jS] 
and Eq.[6l In Table [5] all seasonal fits for differential rotation are 
listed. 

Because most spots on XX Tri appear at high latitudes, Eq.[6| 
with its sin^ b term leads to a much better fit to the observed 
shear than Eq. [S] Therefore, we favor Eq. [6] over Eq.jS] All ob¬ 
serving seasons show a solar-like differential rotation law with 
an overall shear parameter of a = 0.016 + 0.003 and a lap time 
of 1500 days. The surface shear on XX Tri is therefore only 
around a tenth of the solar value. 


Because turbulent diffusion is believed to be the dominating 
effect of spot decay, the decay rate of s pot area is directly 
proportional to the tur bulent diffusivity (iMever et al.l 119741: 
ikrause & Rudigej[l975l) . 

dAjdt = -AnrfT . (7) 

Using the mean decay rate of D = -0.022 + 0.002 SH/day 
from our analysis in Table [4] leads to a turbulent diffusiv¬ 
ity of rjT = (6.3 + 0.5) X 10^"^ cm^s. This value is at minimum 
one order of magnitude higher than that predicted for solar 
value s, which vary f rom 10^^ cm^/s (iDikpati & Charbonneaul 
1 19991) to 10^^ cm^s (iRudiger & Kitchatinovll2000l) . The diffu- 
sion timescale for the magnetic field inside the convection zone 
(CZ) is given by 


rjT 

where the width of the stellar convection zone. Using stel¬ 
lar models calcula ted with the Yale R otational stellar Evolution 
Code (YREC; see lSpada et al.l (|2013|) for more details) we esti¬ 
mate a depth of 0.94 for the convection zone of XX Tri. Thus, 
leading to a magnetic cycle of approximately 26 + 6 years. 

5.5. A starspot movie 

We merged all 36 Doppler images into an animated gif file 
(xxtri-di-anim.gif), which is available through the A&A 
video depository or our group web pag4!l- It particularly empha¬ 
sizes surface detail not covered in the spot-decay analysis and 
shows the general migration and decay trends. It may be used 
for general demonstration purposes. 

The movie shows the same four equidistant surface views 
as, e.g., in Eig. [6]but just as a function of time. Between each 
Doppler image a time delay of 250 ms is included. 


6. Discussion and summary 

Thanks to our robotic STELLA telescopes the present time se¬ 
ries of Doppler images resolves 36 single stellar rotations over 
a time span of six years. The sample is long enough that it en¬ 
ables, for the first time, a direct determination of a starspot de¬ 
cay law. Our target is the KO giant XX Tri with a rotation period 

24 d) comparable to that of the Sun but being more massive 
by 26 % and significantly older with an age of 8 Gyrs. This 
combination of parameters is only possible because the star is a 
component of a close binary. A comparison with solar analogies 
is therefore only for general guidance. 

The time series enabled the cartography of a variety of sur¬ 
face activity phenomena, such as active longitudes, fiip-fiops, 
and differential rotation on XX Tri. However, our main re¬ 
sult is a spot decay law leading to a prediction of a mag¬ 
netic activity cycle solely based on an observationally con¬ 
strained value of the turbulent magnetic diffusivity. The indi¬ 
vidual spot decay rates, dAjdt = D, scattered between -0.019 
and -0.025 SH/day over the six year observing period with a 
mean value of D = -0.022 + 0.002 SH/day. The rates for spot 
growth were between +0.013 and +0.030 SH/day with a mean of 
D = +0.021 + 0.002 SH/day and thus of nearly the same amount 

^ http://www.aip.de/en/research/research-area-cmf/ 
cosmic-magnetic-fields/stellar/stellar-activity/news 
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Table 5: Results for differential rotation. Parameters for best fits using Eq. [5] (middle part) and Eq. [6] (right part). The parameter Qeq 
is fixed for P^ot = 24.0 d in both cases. 


Season 

# ccf maps 

AQ (7d) 

a 

lap time (d) 

1 Ui (7d) 

Uz (7d) 

a 

lap time (d) 

2006/07 

6 

0.05+0.03 

0.003+0.002 

^6870 

0.53+0.04 

-0.68+0.04 

0.010+0.001 

;^2430 

2007/08 

4 

0.30+0.08 

0.020+0.005 

^^1190 

0.83+0.20 

-1.35+0.23 

0.035+0.010 

^690 

2008/09 

6 

0.29+0.06 

0.019+0.004 

^^1240 

0.51+0.16 

-0.95+0.18 

0.030+0.011 

^810 

2009/10 

4 

0.11+0.04 

0.007+0.003 

^3290 

0.22+0.21 

-0.39+0.24 

0.011+0.012 

;^2130 

2010/11 

3 

0.19+0.06 

0.013+0.004 

^1870 

-0.20+0.25 

0.01+0.33 

0.013+0.058 

^1890 

2011/12 

6 

0.00+0.04 

0.000+0.003 

— 

0.64+0.14 

-0.74+0.17 

0.007+0.002 

;^3460 

2006-12 

29 1 

1 0.13+0.04 

0.009+0.003 

;^2740 1 

1 0.45+0.07 

-0.69±0.09 

0.016±0.003 

^1530 


as the decay. In above units, the spot decay on XX Tri wou ld be 
of the order of 10"^ times faster than what iBuinb^ (1 19631) sug¬ 
gested for sunspots. Bumba proposed a mean value for D of 
-4.2 MSH/day. However, the areal size of starspots on XX Tri 
is also 10^ t o 10"^ times larger than the l argest observed sunspots 

10“^ SH: lBaumann & Solankill200^ . Erom this, we depict a 
turbulent diffusivity for XX Tri of 6 x 10^"^ cm^/s, a value be¬ 
tween 10 to 10,000 times larger than current model values for 
the solar convection zone (from the surface layers to the bottom 
of the convection zone). Because the (squared) absolute depth of 
the convection zone of XX Tri is about 1,200 times larger than 
that of the Sun (200 Mm), the diffusion timescale becomes com¬ 
parable to that of the Sun. We obtain an average diffusion time 
of 26 yr for XX Tri compared to 12 yr for the Sun (the latter 
for an assumed diffusivity of 10^^ cm^s). 

So far, stellar activity cycles were inferred from long-term 
chro mospheric Can H&K or photospher ic V-band variations 
fe.g.. lBaliunas et all 1 9951:[Olah et al.ll2007h or from repeated de¬ 

tections of a “fiip-fiop” phen omenon (see lHackman et al.ll2Q13l: 
iKorhonen & Jarvin^ l2007l) . Derived tin iescales for RSCVn 
stars range between 2-50 yr. Just recently, [Olah et al.] (l2014l) in¬ 
vestigated 28 years of photometry of XX Tri along with two 
other over-active K giants and found a long-term sinusoidal 
brightness trend with a length comparable to the length of the 
data set, i.e., 28 yr. Our diffusivity-based cycle prediction 

of 26 yr matches this observation surprisingly well despite 
the comparable shortness of the photometric coverage. XX Tri 
reached its maximum brightness in 2009 an d is in a declining 
state since then. Removing this overall trend. lOMh et al.l (l2014l) 
found a second shorter-period variation of approximately 6 yr or 
an integer multiple of it. No clear explanation for this period, if 
real, could be given. 

As already mentioned in Sect. 1431 the reconstructed giant 
spots may be monolithic, but could also be a conglomerate of 
smaller, unresolved spots. This “classical” uncertainty could in 
principle impact on the interpretation of the observed decay rate, 
and thus the cycle length. Assuming that each unresolved spot 
only decays (and never grows), and does not interact with an 
other spot fragment, then we should observe on average the 
same decay rate as if the spot were monolithic. If decay and 
growth of individual fragments coexist, then our determined de¬ 
cay rate would be just a lower boundary. It then enables a maxi¬ 
mum decay rate of -0.03 SH/day, resulting in a cycle length 
of 19 yr. This cycle period would be close to our determined 
1-cr uncertainty. On the contrary, no such cycle period was de¬ 
tected from broadband photometry. In addition, a possible cycle 
length of six years (from photometry) or two years (from flip- 
flop) would suggest a decay rate of -0.10 and -0.29 SH/day, 
respectively; neither of which is supported by our analysis. We 


conclude that our predicted cycle of 26 yr, including an uncer¬ 
tainty of 6 yr, appears to be the most credible. 

Active longitudes are a common feature in rapidly-rotating 
active stars and there is even some evidence for long-term 
active longitudes and a 7-y r flip-flop period on the Sun 
(iBerdvugina & Usoskinll2003l) . Our Doppler imagery provides 
evidence for a 2-yr flip-flop period on XX Tri with a preferred 
longitude typically facing the (unseen) companion star. This is 
significantly shorter than the 6-yr period from photometric data, 
but could be the true flip-flop cycle length and thus would iden¬ 
tify the 6-yr period just as an alias. Erom a theoretical perspec¬ 
tive, flip-flops possibly represent the nonaxisymmetric compo¬ 
nent of a mixed-mode dynamo for weakly di fferentially rotating 
stars (lElstner & Korhonenll2005l:lMossll2004l) . These authors ex¬ 
plained flip-flops as an excited nonaxisymmetric dynamo mode, 
giving rise to two permanent active longitudes in opposite stellar 
hemispheres, but still need an oscillating axisymmetric magnetic 
field in parallel. The stability of this kind of a mixed mode is still 
a matter of discussion. 

Eor XX Tri, our Doppler images indicate a weak solar-like 
differential rotation of a = 0.016 + 0.003, which seems to be 
a typical value for this kind of rapidly rotating stars. Despite 
the fact that large-scale mean-field dynamo models predict only 
a poo r tracing qual ity for its large sp ots (iKorhonen & Elstn^ 
1201 ll: but see also ICzesla et~^ l2013l) . numerous differential- 
rotation laws were deduced from cross correlations of cool fea¬ 
tures from consecutive Dopple r images (at this point we refer 
to the many references cited in IKorhonen & Elstnerll201l1) . Re¬ 
cently, weak solar-like differential rotation was confirmed, e.g., 
on the K giants I E Hva (iKovari et al.l20l4IWeber & StrassmeieJ 
Il998l) or f And (iKovari et al.ll2012ll2007al) . while weak antiso¬ 
lar differential rotation wa s confirmed on the K-giant cr Gem 
(IKovari et al.ll2015l l2007bl) . Eurthermore, weak antis olar differ¬ 
ential rotation wa s claimed for the K gian ts UZ Lib (lOlah et al.l 
l2003l). HD 31993 (IStrassmeier et al.ll200^ . and possibly HU Vir 
(IStrassmeie^l 19941) . All the latter still in need of further indepen- 
d ent verification. Eor a previous summary on this topic we refer 
to lWeber et al.l (l2005l) . XX Tri fits into the differentially rotating 
giants with approximately ten times weaker surface latitudinal 
shear when compared to the Sun. 
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Fig. A.l: Doppler images of XX Tri for the observing season 
2007/08. Each DI is shown in four spherical projections sepa¬ 
rated by 90°. The rotational shift between consecutive images is 
corrected, i.e., the stellar orientation remains the same from map 
to map and from season to season. The time difference between 
each Doppler image is indicated in units of rotational phase 0. 


Appendix A: Doppler images for the seasons 
2007/08 to 2011/12 

Appendix A.1: Season 2007/08 

In Fig. IA.il five almost consecutive Dopple r ima ges are shown, 
which cover around eight rotations. In Fig. IA.2l the spot-model 
fits of the Doppler images are shown. 
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Fig. A.2: Spot-model fits of the Doppler images in Fig. IA.il Each 
spot is shown with different color/contrast for better visualiza¬ 
tion. 


Appendix A.4: Season 2010/11 

In Fig. IA.7I five almost consecutive Doppler images are shown, 
which cover around nine rotations. In Fig. I A. 8 1 the spot-model 
fits of the Doppler images are shown. 

Appendix A.5: Season 2011/12 

In Fig. lA.9l seven almost consecutive Doppler images are shown, 
which cover around ten rotations. In Fig. lA.lOl the spot-model 
fits of the Doppler images are shown. 


Appendix A.2: Season 2008/09 

In Fig. lA.SI seven almost consecutive Doppler images are shown, 
which cover around nine rotations. In Fig. IA.4I the spot-model 
fits of the Doppler images are shown. 


Appendix B: Line profiles of Doppler images for the 
seasons 2006/07 to 2011/12 

Fig. lB.ll and Fig. lB.2l show the observed and inverted line profiles 
of each Doppler image for all observational seasons. 


Appendix A.3: Season 2009/10 

In Fig. IA.5I five almost consecutive Doppler images are shown, 
which cover around eight rotations. In Fig. lA.bl the spot-model 
fits of the Doppler images are shown. 


Appendix C: Phase coverage of Doppler images for 
the seasons 2006/07 to 2011/12 

Fig. IC.ll shows the phase coverage of each Doppler image for all 
observational seasons. 
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Fig. A.3: Doppler images of XX Tri for the observing season 
2008/09. Otherwise as in Fig. IA.il 


Fig. A.4: Spot-model fits of the Doppler images in Fig. lA.31 Oth¬ 
erwise as in Fig. lA.2l 


Appendix D: CCF-maps and DR-fits for the seasons 
2006/07 to 2011/12 


Appendix E: Longitudinai spot distribution for the 
seasons 2006/07 to 2011/12 

Fig. IE. II shows the mean distribution of the spot area from our 
spot-model fits for each observing season. 


Fig. ID. II shows the ccf maps from each observing season. In 
Fig. ID. 21 the observed differential rotation pattern determined 
from the ccf maps, together with the best fit of the differential 
rotation following Eq.[5]and Eq.[6]are shown. 
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Fig. A.5: Doppler images of XX Tri for the observing season 
2009/10. Otherwise as in Fig. IA.il 
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Fig. A.6: Spot-model fits of the Doppler images in Fig. lA.51 Oth¬ 
erwise as in Fig. lA.2l 
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Fig. A.7: Doppler images of XX Tri for the observing season Fig. A.8: Spot- model fits of the Doppler images in Fig. lA.71 0th- 
2010/11. Otherwise as in Fig. IA.il erwise as in Fig. lA.2l 
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Fig. A.9: Doppler images of XX Tri for the observing season Fig. A. 10: Spot-model fits of the Doppler images in Fig. IA.9I 
2011/12. Otherwise as in Fig. lA.ll Otherwise as in Fig. IA.2I 
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Fig. B.l: Line profiles of Doppler images #1-18. Each figure shows the observed (solid lines) and inverted (dotted lines) line 
profiles for one Doppler image stating their mid times and the respective phases. Rotation advances from bottom to the top. Their 
corresponding RMS-errors are given in Table [S] 
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Fig. B.2: Line profiles of Doppler images #19-36. Otherwise as in Fig. IB. 11 
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Fig. C.l: Phase coverage of Doppler images from 2006 to 2012. Different filled (colored) symbols represents the phases of each 
individual Doppler image, whereas not-filled circles represents non-used spectra (except for gap filling). The arrows indicate the 
spectra, which were used to fill up large observational gaps. Detailed information is given in Table [S] 
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Fig. D.l: Cross-correlation function maps from 2006 to 2012. Each map represents the average ccf map for one observing season. 
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Fig. D.2: Differential rotation signatures from 2006 to 2012. Analyzing the ccf maps in Fig. ID. li re veals a weak solar-like differential 
rotation. The dots are the correlation peaks per 5°-latitude bin and their error bars are defined as the FWHMs of the corresponding 
Gaussians. The dashed line represents a fit using Eq.[5l whereas the solid line represents a fit using Eq.[6l The parameters for each 
fit are summarized in Table [S] 
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Fig. E.l: Longitudinal spot area distribution on XX Tri from 2006 to 2012. Shown are the seasonal mean distributions of the 
individual spots (solid colored lines) from our spot-model fits for each observing season. The black dashed line represents the total 
spotted area. The spot area is given in solar hemispheres on the left axis (1 SH = 3.05 Gm^) and relative to the total area of a stellar 
hemisphere of XX Tri on the right axis. 
















































